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We present the differential rates for exclusive B —>■ X£i£ 2 , where is a charged massless lepton 
and £2 is a charged or neutral massless lepton, and X is a mesonic system up to spin 2. The cases 
of interest are semileptonic (SL) B —>■ X^c,u}^~£'e decays, and B —>■ Xs£~£'^ where the the di-lepton 
can be cc resonances or non-resonant electroweak penguins (EWP). We consider helicity amplitudes 
having non-zero relative phases that can be potential new sources for CP-violation. Our motivations 
for these additional phases include a complex right-handed admixture in the hadronic weak charged 
current for the SL decays and complex Wilson coefficients in the effective Hamiltonians for the EWP 
decays. We demonstrate the efficacy of a novel technique of projecting out the individual angular 
moments in the full rate expression in a model-independent fashion. Our work is geared towards 


ongoing data analyses at BABAR and LHCb. 
PACS numbers: 12.15.-y,12.10.Dm,13.20.-v,12.15.Hh 

I. INTRODUCTION 

The theory of semileptonic (SL) B decays is a rich and 
well-studied subject [1-6]. Within the framework of the 
Standard Model (SM), this has been widely used to probe 
the nature of the electroweak interaction and the struc¬ 
ture of the Cabibbo-Kobayashi-Maskawa (CKM) ma¬ 
trix. In particular, the CKM matrix elements \Vcb\ and 
|I4b| can be extracted from the rates of the processes 
B -A where X^c,u} is an exclusive charm or 

charmless meson state, respectively. 

The full differential rate in the SM for these processes 
have been previously presented by several authors in 
Refs. [1-6]. The current article extends these results in 
the following fashion. Instead of assuming the relevant 
helicity amplitudes to be relatively real, as is the current 
status, we provide expressions corresponding to complex 
amplitudes. A specihc motivation for admitting com¬ 
plex amplitudes in SL decays is to consider, instead of a 
purely left-handed (LH) weak charged current as in the 
SM, an additional complex right-handed (RH) admix¬ 
ture, cr, that could arise in new physics (NP) scenarios, 
as shown in Fig. 1 [7, 8]. A complex non-zero cr leads to 
additional angular terms in the full differential rate. In 




FIG. 1: (Color online) Effective RH coupling arising in (a) 
gluino loop in supersymmetric models and (b) heavy RH W 
boson in left-right symmetric models. 


particular, a non-zero phase in cr can lead to CP viola¬ 
tion in SL B decays [4, 8-10]. 

Consider on the other hand the process B -A Xi~ , 
where X subsequently decays into two pseudoscalars, and 
the charged di-lepton system can be either be a cc res¬ 
onance (J/'i/’i V^(2S')) or non-resonant electroweak pen¬ 
guins (EWP). It is well known that the helicity ampli¬ 
tudes here have non-zero relative phases [11]. Compared 
to the SL case, where the leptonic current is purely LH, 
both LH and RH components exist for the charged di¬ 
lepton case. The LH and RH terms add incoherently to 
give the total rate. Therefore, while the number of an¬ 
gular observables remain the same, the number of inde¬ 
pendent real amplitude components to extract increases 
almost two-fold. The angular observables are not inde¬ 
pendent which leads to ambiguities in the solutions for 
the amplitudes [12, 13]. 

A simplification occurs for the case where the di-lepton 
is a resonant cc meson that decays electromagnetically. 
Since electromagnetic interactions conserve parity, the 
LH and RH amplitudes are equal for this case. The re¬ 
duced number of real amplitude components result in a 
single two-fold ambiguity, as explained in Sec. VIII. 

To sum up, in this article, we examine the generic 
B -A Xixii decay, where Li is a charged massless lep¬ 
ton and £2 is a charged or neutral massless lepton, and 
both the LH and RH helicity amplitudes can be non-zero, 
complex and independent of each other. The SL and the 
resonant cc instances represent special cases leading to 
certain simplifications. We expand the full angular ex¬ 
pression in an orthonormal basis of spherical harmonics 
and provide moments to project out each angular com¬ 
ponent. Since the basis is orthonormal, this reduces to a 
simple counting measurement. We explain how to extract 
the covariance matrix of the moments and the treatment 
of background subtraction, again, as counting measure¬ 
ments. As long as the set of basis functions is “large- 
enough” , our method is the most model-independent way 
of describing the data, as inputs to theory modeling. 
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FIG. 2: (Color online) The quark level diagram for the SL FIG. 3: (Color online) The set of four kinematic vari 
decay B(bq') ^ X{qq')£-Ve in the SM. <(> € {g , x, cos 0^, cos 0v} for the SL decay chain 

i? —> p(—> 7r7r)tT*(—>■ £ ve)- 


II. THE KINEMATIC VARIABLES 

Consider the SL decay process B{bq') —)■ X{qq')£~V^ 
shown in Fig. 2. At the quark level, in the SM, this is 
a flavor changing process where a heavy b quark emits a 
charged W* (off-shell) and decays into a lighter q G {u, c} 
quark, with the decay vertex containing the CKM matrix 
element Vqb- An important feature of SL B decays is that 
the leptonic side interaction is well-understood, thereby 
facilitating study of the complicated non-perturbative 
QCD interactions that reside on the hadronic side. The 
momentum transfer squared between the leptonic and 
hadronic systems is The hadronic side is thus probed 
by the dependent form-factors (FF), just like in deep 
inelastic scattering (DIS), save that g^ > 0 is now time¬ 
like, instead of spacelike in the DIS case. For the EWP 
case, the W* can effectively be thought as being replaced 

by{ 7 *,^*}- 


A. Kinematics 

Without loss of generality, we take £i = l~ and 
I 2 = Vi- We denote the 4-momenta of the parent B, 
the daughter meson A, the charged lepton £~ and W as 
Pb, Px, Pe and respectively. The W* 4-momentum is 
q = {pb -px), so that 

q^ = {PB - Pxf' = m%+mx - imsEx (la) 

Pb Px Ex 

W = Vb ■ Vx = - ■ - = - 

ms mx mx 

= + (lb, 

2mBmx 

where Ex is the energy and w is the 7 -factor of the X as 
seen in the B rest frame (RF). If we consider the breakup 
of i? —>■ XW* as a two-body decay, where the virtual W* 
boson has mass the two-body breakup momentum 

is given by 


virtual mass, ^/^ = \/q^^ = {ruB — mx)- Since the 
meson X is at rest in the B RF now, the 7 -factor w = 
Wmin = 1- This kinematic region is convenient for lattice 
and heavy quark effective theory calculations. On the 
other hand, at g^ = g^j^^ « 0 (for the massless leptons), 
the breakup momentum is largest 


k 


max 


9 9 

TUq — 

2mB 


( 3 ) 


Since the breakup momentum and the 7 -factor w are 
related as 


k = mx — 1. 


we also have 


771 ^ + mx 

2mBmx 


or, the /3-factor as 


P 


max 


-m\ 
mpg -I- m^x 


( 4 ) 

( 5 ) 

( 6 ) 


This “large-recoil” region is convenient for light-cone sum 
rules and soft collinear effective theory calculations. 

When the outgoing meson in a vector meson, its po¬ 
larization is important as well. The vector meson decay 
products act as the analyzer. For example, in the case 
of p —>■ TTTT shown in Fig. 3, the analyzer (A) is the tt 
momentum direction in the p helicity frame with respect 
to the B RF. This defines the helicity angle Oy- For 
uj —> 7 r+ 7 r“ 7 r*^, the normal to the w decay plane plays the 
role of the analyzer. The last additional kinematic vari¬ 
able is Xj tbe dihedral angle between the W* and the 
vector meson decay planes in the mother B RF and care 
must be taken to note the quadrant of the angle x (see 
Fig. 4). We refer to the set of four kinematic variables as 
(/) = {g^, cos cos 6 *y,x}- 


Two kinematic limits are of special interest. At “zero- 
recoil” , k = 0 and the IT* attains the maximum allowed 


B. Sign conventions of {0i,0v,x} for B Xiii 2 

We stress here that the angles in this section are for the 
B decay (that is, a b quark transition). The CP conjugate 
case is decsribed in Sec. VI. 
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FIG. 4: (Color online) The reference frames for calculating 
6t, 6v and x for B —>■ VW*{-^ £~i'e) conforming to (a) Fig. 2 
in Gilman-Singleton [1] and (b) Fig. 3 in Hagiwara [5]. The 
relations between the two sets of angles are given in the text. 


1. SL case: {£ 1 , 12 } = {£ 


We follow the definition of the angles in Fig. 2 of 
Gilman-Singleton [1]. We first boost everything to the 
B RF. There are two sets of co-ordinate axes, {xg, yi, Zi} 
and {xv, yv, w}, as shown in Fig. 4a for the vector me¬ 
son (V) case. These are the helicity frames of the W* and 
the V. The connection is that xi = —xy, yi = yv and 
zi = —Zv- The dihedral angle X = V’e + V’V: where we 
note that the azimuthal angles :pi and ipv are calculated 
in two different frames. We set = 0 by ensuring that 
the charged lepton £ lies in the x^-Zi plane and and has 
the a:-component of its momentum > 0. This completely 
fixes the quadrant of y, and therefore the signs of its sine 
and cosine. 

To measure 9( and Oy, we boost to the W* and V 
rest frames and measure the polar angles of the £ and 
A, respectively. Here A is the analyzing direction of the 
vector meson decay as tabulated in Table L 

Korner-Schuler [3] and Hagiwara [5, 6 ] follow a differ¬ 
ent convention where both the orientations of the axes 
for the leptonic and hadronic systems are the same. The 


Decay topology 

A direction 

D* - 

> Dn 

Pd 

D* —> D 7 

Pd 

P'^/fo - 


PnO 

P 7/0 -s- tt+tt 

Pn + 

u) — >■ 7r"'"7r“7r‘^ 

Pn+ X Pn- 

K* Kn 

Pk 


TABLE I: The analyzing direction A in Fig. 4 for the different 
mesons in a B decay. See also Sec. VI for the B case. 


relations are 

0KS ^ (7^) 

^ (7b) 

^KS ^ ^ + ^GS (7^) 

where the “KS” superscript refers to Korner- 
Schuler/Hagiwara and the “GS” superscript to 
Gilman-Singleton. 

The conventions followed by Richman-Burchat [2] 
(“RB”) on the other hand are related to the GS defi¬ 
nitions as 



Ill 

(8a) 

iRB 

V 

III 

Ol 

(8b) 

,RB 


(8c) 


We adhere to the GS conventions in this work. 


2. EWP case: {^,£ 2 } = {£~ 

We again follow the GS conventions with the sin¬ 
gle replacement vi —>■ everywhere. Compared to 

other EWP conventions in the literature [13-18], the only 
change is 

where the superscript “EWP” refers to the aforemen¬ 
tioned theory references (see also appendix). 

III. EFFECTIVE HAMILTONIANS 
A. SL decays 

Consider the process b —> q£~ve (where q € {c,u} and 
£ € {e,/r}) in terms of an effective 4-Fermi interaction 
Hamiltonian: 

+q { 9 ^ + 9^lb) hivh] , ( 10 ) 
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where we have assumed only LH neutrinos {i'l = 
and neglected any tensor terms associated with baryon 
and lepton number violations (leptoquark models [19]). 
Here, = Vqb denotes the usual LH CKM matrix el¬ 
ement in the SM. The vector and the axial interactions 
are written as gy = {^ + en) and gA = [^ — £r) and 
in general, cr is allowed to be complex to incorporate 
additional CP violating effects [4, 8-10]. There are also 
two terms, g^ and g^, corresponding to scalar and pseu¬ 
doscalar interactions, respectively. To retrieve the SM 
part, one puts cr = g^ = g^ = 0 . 


1. The B ^ SI Vi case 

The transition matrix element pertaining to the pro¬ 
cess B —>• S£~Vi, where S' is a a 0+ scalar state, is 
then Ms = (SlTLeslB). We note here that a nega¬ 
tively charged lepton and a RH anti-neutrino is being 
produced (since we have allowed for extra phases, we have 
to be careful about CP conjugates now). The hadronic 
matrix elements corresponding to the terms gM^,TS} 
are written in terms of two form factors u+{q'^) and 
Uo{q^) [20, 21]: 

{S\Q7fj.^\B)v = 0 ( 11 a) 

{S\qJt,j5b\B)A = u+{q'^)(^{pB +ps)v - 

, i2^iPB+Ps)■q 
+ uo{q )-^( 11 b 

{S\qb\B)s = 0 ( 11 c) 

(S|g75b|j3)p » Uo{q^) ^^ . (lid) 

mb + mq 

Since parity factors multiply, the right hand side in 
Eq. 11a has to be an axial vector, which one can not 
construct out of the two vectors ps and pb- Therefore 
only the g^ term survives in Eq. 11b, while the g^ term 
in Eq. 11a is zero. Eq. 11b has been written in a form 
that is non-singular at q^ —> 0. However, for the light 
leptons, the q^ terms in go to zero when dotted with 
the leptonic charged current This can be seen by 

using q = (pi+Pv) and the Dirac equation for the (mass¬ 
less) leptons. Hence, all terms proprtional to q^ can be 
dropped. Eq. 11c and Eq. lid follow from Eq. 11b and 
Eq. 11b, respectively, by dotting with q^ and invoking the 
Dirac equation at the quark level. In all, the transition 
matrix element reads: 

Ms = \gAU+{q^){pB+Ps)iih^vL 

+ ( 12 ) 
nib + mq J 

As we will see later, the gp term can be ignored for the 
massless lepton case, and only the gA term will remain. 


2. The B ^ Pi Vi case 

When the outgoing meson is a 0“ pseudoscalar state P, 
following the argument given above, the gA term vanishes 
and there is only a gv contribution, with the two form 
factors /+(g^) and /o(g^): 

(F|g7^6|i?)y = /+(g^) ^(pe -bpp)^ - ^ 

q 

(13a) 

{P\qTM5b\B)A = 0 (13b) 

(P|g 756 |H)p = 0 (13c) 

{P\qb\B)s^foiq^ r"\'^"^ , (13d) 

nib + mq 

and the amplitude reads: 

Mp = {gv/+(g^)(PP + Pp)bih^BL 

+ (14) 

mb + mq ) 

As in the scalar case, the gs can be ignored for the mass¬ 
less lepton case, and only the gy term will remain. We 
note that the structure of of Eq. 12 and Eq. 14 are quite 
similar, except for the coupling terms and the form fac¬ 
tors. Since gA and gy are proportional to (1 — ep) and 
(l-bep), respectively, the effect of a non-zero ep is differ¬ 
ent between the outgoing scalar and pseudoscalar meson 
states. 



3. The B ^ Vi Vi case 


When the outgoing meson is 1“ vector meson E, both 
the gy and gA terms contribute and the hadronic current 
can be written in terms of four form factors Aq, Ai, A 2 
and V: 


(Vjqj^bjBjy = (15a) 

my + mp 

{y\qiiil^b\B)A = 2TOyAo(g^)^^y^gp 
-b (mp -b my)Ai{q^) 


{V\q-fbb\B)p 


A2{q^ 




mp + my 

{pb +Pv) 


PB +Pv)i 


2myAo{q‘^) 


q" 

£y • q 

mb -b mq 


-q^i 


{V\qh\B)s=0, 


(15b) 

(15c) 

(15d) 
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and the matrix element is: 


2GFVgbf\ f 2iV{q^) ,, 

■^V — - 7 =- ^ 9v I - 1 - ^iivaP^V PvPb 

V2 I L \'^v +mB 

- 9^(^{mB + mv)Ai{q^)ey^- 

h^VL 

IV ' _ 

2m B 


^2(g^) [PB +Pv)^ 

mB + Mv 


mb + mq 


Ey ■ q Ao{q )] evLf ■ (16) 


B. EWP decays 

The effective Hamiltonian for 5 —)■ s transitions can be 
expanded in the form [14, 15, 22] 


tree-level operators, while 03 <i<Q are gluonic penguins. 
Os is a gluonic dipole operator and is also suppressed by 
a power of The scalar and pseudoscalar operators. 

Os,Pi do not contribute in the SM. 

The three main contributing terms for the b —)■ s£~£'^ 
EWP decays are O 7 , Og and Oig. O’j is the 7 * penguin, 
while Og and Oig get contributions from the Z* penguin 
and W* box-diagram. 

The Wilson coefficients are calculated by matching the 
effective and full theory at the /i ~ mw scale, and evolved 
down to /i ~ mb by the renormalization group equations. 
In the SM, the rough hierarchy is C 7 ~ —0.331, Cg ~ 
4.27 ad Cio ~ —4.173, so that Cg and Cig contributions 
dominate, except at q^ —> 0 , where the photonic penguin 
dominates. 

Following Ref. [18], we next define the following coef¬ 
ficients 


E E (17) 

E q ^ 

where q G {u, c, t} represents the quarks running in 
the loop (dominated by the heavy top quark) and i G 
{1,..., 10, S', P}. The unprimed and primed components 
represent the LH and RH (absent in the SM) hadronic 
currents, respectively. The C^’s are the scale depen¬ 
dent Wilson coefficients that encode the short-distance 
physics, while the Oi’s are local operators representing 
the non-perturbative long-distance physics. The explicit 
forms of the operators can be found in Ref. [22]. Oi<:e are 
the 4-quark operators, suppressed at leading-order in the 
SM, but can contribute via charm-loop effects, especially 
near the charm-threshold in q'^. Out of these, Oi <:2 are 






1^9 ^ 9/ + \^10 ^ 10/ 

(Cf+ C«%)T(Cio + C<^®'io)l 


/2 

/2 


C 7 = [Cf - C^%)/2 
C\ = {Cf + C''=®7)/2, 


(18a) 

(18b) 

(18c) 

(18d) 


where {P, P} signifies the handedness on the leptonic side 

and the expressions of C®** can be found in Ref. [15]. 
It should be noted that the effects of charm loops (from 
Ci<g) enter Cf, thereby incorporating strong phases 
into 

For X being in a spin-J state, the helicity amplitudes 
in terms of the < 7 ^ dependent form-factors Fi,t, ^ 1 , 2 , V 
and ri ,2 are [16]: 


H 


L,R 


2 mBk j-^L,R 
J=0~ { 


i C^’“Fi(g2) + ^7 Ffq^) 


rriB + mx 


H. 


L,R 




= /3j , 

j>i \mx 


C^’^^mB + mx)Ai{q‘^) + ^^(m| - m\)C7T2{q^) 


(19a) 


=F 2 TOBk 


tnB + mx q 


H, 


L,R 


aj 


( — ) 
J>i 2mx\f(f \iTixJ 


j-i 


^l,r 


( to | - mx - q^){mB + mx)Ai{q^) -—^ 2 ( 9 ^) 


2TnBC7 


(m| -h 3 to| - q^)T 2 {q^) - 


ruB + mx 

dm^k 
rug — m^^ 


-Tsiq^) 


(19b) 


(19c) 


The {aj,/3j} factors come from Clebsch-Gordon coeffi¬ 
cients and are {1,1} and {y^2/3, l/-\/2} for the vector 
and tensor states, respectively [16]. The k'^”^ terms are 
additional kinematic factors incorporating the angular 


momentum barrier factors for higher spins (see also dis- 
cusion in Sec. V). 

Note that in the above equations for the helicity ampli¬ 
tudes, relative to the convention in the EWP literature, 
we have taken out an overall normalization factor. The 
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benefit is that the SL limit is easily arrived at by the sub- 
stititions = C'^ = Cr = C'-j = 0, and = I 

(see also Eq. 23). The terms corresponding to Ft and 
Ti's do not exist in the SL case and Fi is identified as 
M+(/+) for B S{P)£~Pe- 


For the scalar term L 5 , the helicities of the two leptons 
must be the same, since the total spin of the di-lepton 
system is 0. This means, that although [5, 6 ] £vl = \/^, 
so that 

Ls = 2lvL = 2y^, (26) 


IV. DIFFERENTIAL RATE FOR B Xiih 
Following Hagiwara [5, 6 ], the differential rate is 

dr=^ E (20) 

2 tob 

nnai spins 

where the incoherent sum is over the spins of all final- 
state particles and three-body X£i £2 phase-space factor 
is 


k dq'^d cos 9^ 
= 2mB (47r)3 


( 21 ) 


where k is the usual X 3-momentum magnitude in the 
B RF. Including only spin-0 and spin-1 states for the di¬ 
lepton system, the invariant amplitude can be written in 
the form 


\ + '^s)Ls + XI 51 \ > 

^ [ r/=±l Ae{o.±i} J 

( 22 ) 

where A is the helicity of the hadronic system X and 
the handedness 77 = (A^^ — A^^) = —1(-|-1) for LH(RH) 
leptonic currents. For SL decays with f 1 always being the 
charged lepton, for a {£~Vi) final state, since A = -1-1/2 
for the purely RH Vi, we have 77 = —1. For the v) 
final state, we have 77 = -1-1 for the purely LH v. Here, 

F is a scale factor that equals GpVqb in SL & —)• <7 type 

_ _ 

transitions. For the EWP decays B —> K £“£+, the 
effective replacement is [18]: 




EWP 


(23) 


The hadronic helicity amplitudes for B (that is, con¬ 
taining a b quark) are defined as [4] 

nl = {e*^,{X))^{XiX)\J^\B)^, (24) 

with the spin-quantization axis along the X flight direc¬ 
tion in the B RF, while the leptonic helicity amplitudes 
are (for massless leptons) 

Ll = 2 (£w* (-A))^ = 2^dl^ie,) (25) 

where the spin-quantization axis is along the di-lepton 
flight direction in the B RF, opposite to the X direction. 
Since the parent meson is spin-0, the helicities of the 
daughter hadronic and the leptonic systems have to be 
the same. 


the lepton spin-configurations for the L\ and Ls cases 
are different and the two components must add incoher¬ 
ently in total rate expression. The spin-0 leptonic current 
terms for the massless lepton case are therefore second 
order corrections relative to the SM, and will be neglected 
henceforth. 


A. The B —>■ P{S)i vi case 

Following the calculations in Ref. [5], one can show that 
the amplitude in Eq. 14 for the SL outgoing pseudoscalar 
meson case is 

Mp = 9v (-4777,Bk sin 9if+{q^)), (27) 

where we have neglected the gs term because, as ex¬ 
plained above, it is a small second order correction to 
the SM. Plugging this into Eqs. 20 and 21, we get 

^ = (28) 
where dcj) = dq^dcosOi. 

The outgoing scalar meson case is obtained by replac¬ 
ing gy with gAi and /+ with t 6 _|_. 


B. The X Pi P 2 case 

For the case where the X system decays into a pair of 
spin -0 pseudo-scalars the amplitude can be expanded as 


j : iMuj 

[r)=±l Ae{0.±l} J 

where, for the hadronic system, 

HlY,^ = Y,Hl’^Yj\9y,x). 

J 


(29) 


(30) 


The spherical harmonics are given in terms of the Wigner 
d-matrices as 

Y^iOy, x) ^ \j‘^-^di,{ey)F^^ (31) 

and the differential phase-space element is now 
d<f> = d(j>3d cos 9ydx 
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Putting everything together, the full expression for the the same form as /if, except with all the LH amplitudes 

4-D differential rate for X decaying to two pseudoscalars replaced by their RH counterparts, 

is then 


dr _ ^ 

d<j) m\{d'KY 


(32) 


C. The B —l VI vi case 


where. 


\M\^= E 

77—±1 


Ae{o,±i} J 


(33) 


and jg relevant branching fraction (BF). 

The LH and RH contributions add incoherently since the 
final-state spin configurations are different on the lep- 
tonic side. 

For J e {0,1, 2}, we denote the spin-0, spin-1 and spin- 
2 helicity amplitudes as re¬ 

spectively, where the superscripts denote the handedness 
of the leptonic current. The full expansion of | Alp yields 
41 angular terms: 

1^1' = 1 ^ E(^^ + (34) 

i—1 

as tabulated in Table H. Here, = ±1 is a sign 

factor dictated by the behavior of the angular part under 
9i ^ TT + Oi, since d\^^{9t,) = d\ _^{'K + 9i). hf is of 


The V P 1 P 2 case is the same as Eq. 32, with only 
the "Hf amplitudes contributing. For V —)■ P 7 , we need 
to incoherently sum over the outgoing photon helicity 
A-,, = ±1 cases separately: 


dr _ 

d(j) m^(47r)4 

A.^=±l Ae{0,±l} 

^ 3GF\Vgb\‘^kq^B^^P'^ 

32m^(47r)‘* 

X E |2A.ySindy'Hf (—sind^) 
A-y—itl 


+ 77 ^ 1(1 -I- A.^ cos6iy)(l — cos9i)e^^ 

+ Tlf 1(1 — A.y cosdy)(l -I- cos9g)e~^^ 

(35) 

where the extra factor of ^ ensures normalization to the 
appropriate BFs. 

We write "Hf = and set do = 0. For B —> 

V£~Vi, the expressions in Eqs. 32 and 35 can then be 
summarized as: 


dF 

d(j) 


C 

l + {l-a )/2 


I idf (1 — cosd^)^-|-idp(l-|-cosdf)^ (1 — a COS^ dy) -f 4idg sin^ a COS^ 


-I- 2aido sin 9t sin 2 dy 


H+{1 — cos9i) cos(x -|- d+) — id_(l -I- cosdf) cos(x — d_) 


— 2 aid 4 .id_ sin^ d^ sin^ dy cos( 2 x + (d+ — d_)) 


(36) 


where a is -1 for V —■ P7 (such as D* —> or a; — >■ 7r7) 
and -1-1 for V —>■ P 1 P 2 (such as p —>■ tttt or D* —)• Dn) 
type decays. The pre-factor term is 


C' = 






mr 


(37) 


For the V P 1 P 2 type cases, Eq. 36 above agrees with 
Eq. 2.20 in Ref. [1]. It also agrees with Eq. 113 in Ref. [2] 
after taking into account the change in the x definition 
as given by Eq. 8 c. 


where the term B accomodates any BEs from the vector 
meson decay chain onwards. 
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1 

6(|77^p + + 8|5'^p + 10|D^p - 8V5Re{S^Dt) 

~\~{L R) 

2 

cosev[12V5Re{H^D^* + + lQVSRe{S^Ht) - Re{DliHt)] 

55 

3 

cos^ ev[S0{\D^\‘^ + |T»£p) - 6(|f7^p + ILffp) + 24|7/o^p + 24^5 Re{S^Dt) - 60|L»fp] 

55 

4 

cos^ 6lv^[-12V5 7?e(7/^H:5:* + 24^15 i?e(i7oE»o*)] 

55 

5 

cos'* ev[-30{\D^\^ + \D^\^) + 90|E»^p)] 

55 

6 

cos^ 6»46(|77^p + - 8|5'*'p - 10|L»^p + 8^7?e(S'^D^*)] 

55 

7 

cos*‘6»<:cos6iv[12^i7e(77^D^:* + - 16^7?e(S'i7o*) + 8^7?e(Do77o)l 

55 

8 

cos^e^cos** 6 Iv[ 30(|D^P + |L)£|2) -6(|77^p + |i7f p) - 24|77o^|2 - 24V5 7?e(S^D^*) + 60|L»f P] 

55 

9 

cos^^ 6»<. cos® ev[-l2V5 ReiH’i^D^* + - 24^15 Re{H^D^*)] 

55 

10 

cos®6»<:cos* ev[-80{\D^\^ + |D^|®) - 90|D^|®)] 

55 

11 

sin6»^cos6i< cosxsmev[-8VSReiiH^ + H^)S*) + 4V^Re((H^ + Ht)D!^*)] 

55 

12 

sin6»<.cos6»<;cosxsin6»v cos ev[-24 Re{{H^ + - 8^^ Re^D^i. + D^)S^*) + 20^3 ReUD^^ + D^)D^*)] 

55 

13 

sin Be cos Be cos x sin Ov cos® Ov [-24^ Re{{D^ + )H^* - 12^15 Re{{Hl + H^)D^*)] 

55 

14 

sin0^ cos0^ cosxsinflv cos® 6'v'[—60\/3 7?e((Il+ + D^)Dq*)] 

55 

15 

sin Oe cos Oi sin x sin Sy [8^ /m((77^ - Hi;)S*-*) - 4^15 7m((H^ - Hi:)Hi'*)] 

55 

16 

sindi cosOi sinxsinev cosev[24Im{{H^ - Hi')Ho*"‘) + 8^15 7m((7)^ - Dt)S^*) - 20V3 7m((7)^ - 7)i:)Hi'*)] 

55 

17 

sin0^ cos0^ sinxsinSy cos® 6v[24\/5 7m((71+ — D^)Hq*) + 12vT5 7m((H+ — H^)Dq*)\ 

55 

18 

sin0^ cos0^ sinxsinSy cos® 6y[60\/3 7m((7)+ — D^)Dq*)\ 

55 

19 

sin® Oe cos 2x[-12 Re{HlH^*)] 

55 

20 

sin® Ot cos 2x cos Ov [-12^5 Re{H^D^* + 77^Hi;*)] 

55 

21 

sin® Be cos 2xcos® 6iv[-60 7?e(77^7)i*) + 12 He(H:{:Hi*)] 

55 

22 

sin® Be cos 2x cos® By [12^5 Re{HlD^* + 77:5:Hi*)] 

55 

23 

sin® Be cos 2x cos'* Bv[B>Q 7?e(77+77i*)] 

55 

24 

sin® Bi sin2x[12 Im{H^H^*)] 

55 

25 

sin® Be sin 2x cos By [12^5 Im{H^D^* + DIH^*)] 

55 

26 

sin® sin 2x cos® 6»v[60 7m(77^77i*) - 12 Im{H^Hb.*)] 

55 

27 

sin® Be sin 2x cos® By [-12^ lm{HlD^* + 77^Hi*)] 

55 

28 

sin® Be sin 2x cos* 0^]—60 Im{D^D^*)] 

55 

29 

cos0,[12(|Hi|®-|H*:|®)] 

-(L R) 

30 

cose^ cos6iv'[24V5 7?e(Hi77i* - H:|:77:5:*)] 

55 

31 

cose^cos® By[60{\Dtf - |77^|®) - 12(|Hi|® - |H^|®)] 

55 

32 

cos Be cos® By [-24^ 7?e(Hi77i* - H^77^*)] 

55 

33 

cos61^ cos* 6iv[-60(|77i|® - |77^|®)] 

55 

34 

sine^ cosxsin6lv[8V3 7?e((H^ - Hi)^*) - 4Vl5 7?e((H^ - Hi)77^*)] 

55 

35 

sin6l^cosxsin6lv cos6iv'[24He((H^ - Ht)Ht) + 8^7?e((77^ - Dt)S^*) - 20^7?e((77^ - 77i)77i*)] 

55 

36 

sin6»^cosxsin6»v cos® 6»v[24^7?e((77^: - D^)Ht) + 12^7?e((H^ - Hi)77i*)] 

55 

37 

sin0^ cosxsin^v cos® 0v[6O\/3 7?e((77+ — 77i)77o *)] 

55 

38 

sin6»<. sinxsin6lv[-8V3 7m((H:[: + Hi)S'*'*) + 4^15 7m((H^ + Hi)77^*)] 

55 

39 

sin6i^sinxsin6»v cos6»v[-247m((H^ + Hi)Ho*'*) - 8^7m((77^ + 77i)5'*'*) + 20^37m((77^: + 77i)77i*)] 

55 

40 

sin61^ sinxsinev cos® 6lv[-24^7m((77^ + 77i)Ho*'*) - 12^7m((H:j: + Hi)77i*)] 

55 

41 

sinS^ sinxsin6'v' cos® 60\/3 7m((77:^ + 77i)77o *)] 

55 


TABLE II: The helicity-basis expansion of \A4\^ in Eq. 34. 



9 


V. INCORPORATING MASS-DEPENDENCES 


When the variation of the invariant mass of the X 
system, m = mx, is no longer negligible, Eq. 32 can be 
extended as 


dr 

d(j)dm 



(m) 

m^(47r)4 




(38) 


where p is the mass-dependent breakup momentum of 
X —)• P 1 P 2 in the X rest frame, po is the value of p 
computed at the (dominant) pole mass toq. The overall 
factor of p/po comes from phase-space. 

The amplitudes incorporate a mass-dependent rel¬ 
ativistic Breit-Wigner (rBW) part 


'^'rBW (™) 


BJip,R) _ 1 _ 

\PoJ pj {pq, R) ml - m? - ’ 

(39) 


(J, A; 1, —A|Lb, 0) as: 


s =s^ 

(43a) 

H± = -^H^±^H^ + -^H° 

V2 

(43b) 

""" v!"’ - 

(43c) 


(43d) 


(43e) 

The superscripts on the rhs denote the Lg values, and 
the amplitudes represent the spin-Ps component of the 
corresponding helicity amplitude. Each spin-Pe compo¬ 
nent of the helicity amplitudes acquires a nominal barrier 
factor that scales as b = k^®P^®(k, P). We define the 


normalized quantity xlb = d/ 60 , where we choose to 
calculate the denominator at the pole mass. The mass- 
dependent helicity amplitudes are 


where B'^{p,R) is the phenomenological Blatt-Weisskopf 
barrier factor with R « 0(5 GeV“^), corresponding to 
a meson radius of 0(1 fm). For a P-wave decay this is 
given by [24] 


B'’=\p,R) 


1 


(40) 


For the spin-J resonance having a single decay mode to 
the final state P 1 P 2 , 


r 


J 

total 


= ro 


/ p/m \ 
\Po/moJ 


f B'^(p,R) 
\poJ B-'(po,R) 


(41) 


However, if the spin-J resonance has k decay modes, all 
the individual mass-dependent widths contribute as 

k 

(42) 

where B/ is the branching fraction into the z**' mode. 
Examples of such instances are the decay modes of the 
(/)(1020) or the A|(1430). 

The second form of mass-dependence comes from the 
barrier factor associated with the B decay itself. Let 
the B decay into the di-lepton and X system occur with 
an angular momentum Lg and the break-up momentum 
is k, as given by Eq. 2. If the X system is in spin-J, 
the selection rule is Lg & {J—1,J, J-l-1}. The helic¬ 
ity amplitudes R/ can be re-written in terms of specific 
Lg components with the relevant Clebsch-Gordon factor 


S{m) = xiS 

P±(m)^y ('po 


(44a) 


Ho{m) = 


2 x 0 


± Xi 


Ho + 


H++H. 

2 

H+ - P_ 
2 

H+ + H. 

2 


P±(m) = y (^/3Po + iD+ + P_)) 


H++H_- Ho 


(44b) 

H++H_- Ho 
3 

(44c) 


+ ^iVTb VTb(^+ + ^-)-1/5^0 


± X 2 


10 \^V 10 
D+ - £»_ 


(44d) 


Doim) = ^ (^/3Po + iD+ + P_)) 


■ xi 


'-(Ph.+P_)-. -P 


(44e) 


The mass-independent forms are obtained by the sub¬ 
stitutions xlb 1. For the SL and EWP cases, the 
k-dependent barrier factors were already incorporated in 
Eqs. 19. For the cc —>■ decays, the mass-dependent 

forms in Eq. 44 are more appropriate than the bare am¬ 
plitudes. 


VI. CP CONJUGATION 

Consider the CP conjugation of the process B —)• A(—)• 
PiF)^i^ 2 , where Pi is a charged pseudoscalar meson. 
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and l\ is a charged lepton. The CP conjugate process 
is B > X(—> PlY)^^^ 2 ■ We perform the CP conjugation 
explicitly. That is, for the construction of the angular 
variables, going from B to B, we replace the 4-momenta 
as px -> Px, PPi -t PPi, -> and pt^ p^^. This 
construction leads to x —)■ —x- 

On the other hand, the effect of CP conjugation on the 
helicity amplitudes flips the helicities and weak phases 

'HI{Sw,Ss) ='H_1{-6w,Ss), (45) 

where Swi^s) is any weak(strong) phase and flipping the 
sign of r] changes the LH amplitudes to the RH ampli¬ 
tudes. In the absense of direct CP violation, the simulta¬ 
neous effect of these two transformations is to leave |A4p 
unchanged in Eq. 33. Therefore, with explicit CP conju¬ 
gation of the particles during construction of the angular 
variables (see appendix for details), no additional changes 
to the rate equation are required. 

We also stress here that our unbarred amplitudes are 
dehned for the B (or b quark) decay, in contrast to 
conventions in cc analyses [11, 23], where the unbarred 
(barred) amplitudes are defined for the B (B) decay. 


VIII. THE TWO-FOLD AMBIGUITY 

As mentioned in the introduction, the full differential 
rate does not uniquely determine the helicity amplitudes. 
The ambiguities in the solutions arise from the informtion 
loss in summing over the final lepton spins. A detailed 
study of these ambiguities is beyond the scope of this 
work. However, we point out one particular case. 

Using the identities = (—l)^d:^;^g and d\^ = 
— for ?7 = ±1 and A € {0, ±1}, the expres¬ 

sion in Eq. 33 is seen to be invariant under the following 
global transformation: 

^ (uZl'^y ■ (50) 

We note here again that ry = -|-1(—1) denotes the RH(LH) 
component on the leptonic side. For the electromagnetic 
cc —)■ £~^£~ decays, the LH and RH amplitudes are equal 
and Eq. 50 represents the two-fold ambiguity [11] in the 
determination of /? and /ds from B —)■ J/^/jAT* and Bg —>■ 
J/f/'t/), respectively. 


VII. EXPANSION IN AN ORTHONORMAL 
BASIS 


IX. ANALYSIS FORMALISM 


Equation 32 can be expanded in an orthonormal basis 
of angular functions /i(H) as follows 

^=cx||:/.( 0 )r.(,^)| ( 46 .) 

r,((?2) = rf (g2) + rf (g2), (46b) 

where dO = dcosd^dcosdydx and the super¬ 

scripts in Eq. 46b specify the LH or RH nature of the 
leptonic current. The sign = ±1 depends on the 

signature of fi under Qg —)■ Tr-l-d^. Orthonormality of the 
/i’s imply 

J Mn)fj{n)dn = Sij. ( 47 ) 

The orthonomal angular basis is constructed out of the 
the spherical harmonics = Y™{9£,x) and the re¬ 
duced spherical harmonics P™ = y'^Y^^ {Oy, 0) ■ The 
pre-factor is 

^ =-3m|(4^)4- 

Defining the transversity basis amplitudes j_j as 

ni = {n(^±Hi)/V2, ( 49 ) 

Tables HI and IV list the 41 moments in the helicity and 
transversity bases, respectively. 

We note that since the RH and LH amplitudes are 

equal for the cc —>■ £~£'^ type decays, the terms with 

vanish and only 28 non-zero moments sur¬ 
vive in Tables HI and IV for these cases. 


A. No background case 

1, Method of Moments (MOM) 


Assume a generic rate function constructed out of a 
set of orthonormal basis functions /i(H): 

^^g(0) = ^&J,(H), (51) 

i 

where the aim is to determine the moments bi. We define 
a detector efficiency function e(H), and the normalization 
integrals 


,n) = J (tin) [/,(H)/j(H) • • • fnin)] dn 

'atMC 

ftink)f]ink)■ ■ ■ fnink) , 

fc=i 

(52) 

that are calculated numerically with Monte Carlo 
(MC) events generated flat in dn, and accepted 

events that survive after the detector efficiency is taken 
into account. Also, $ = f dn is the total phase-space 
element. 

The measured moments from the data are 


$ 

JVMC 

gen 


\T d ELX 3. 

2 ’ /* /-7 AT 

h= = I Mn)e{n)—dn = E,,b, ( 53 ) 

k—1 
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i 


ri.hel(^2) 


1 

PSYo° 

[\H^\^ + \Hi\^ + |p£;|2 + |S-^|2 + \d2\^ + \d2\^ + |p£|2] 

+ (_L ^ _R) 

2 

Pi%° 

2 

^ReiH^Bt) + Re{S^H2*) + 1 ReiH^Bl* + HbBb*) 


V 

3 

pOyO 
-*2 ^0 

# + \Dtn - {\H2\^ + \H^n + ^ \H2\^ + ^ + 2 ReiS^Dt) 

V 

4 

pOyO 

[- ReiH^D^.* + H2 d2*) + VZReiH^D^*)] 

V 

5 

pOyO 

1 [-2(|Pi|2 + |P^|2) + 3|Po^|2] 

V 

6 


^ [(|Pi|2 + |P£|2) + + |P^|2) _ 2|5^|2 - 2 |Po^|2 - 2\H2\^] 

V 

7 

P?Y2° 


^ ReiH^B^* + HbBt*) - ^ Re{S^H2*) - | ReiH^Bt) 


V 

8 

P^Y§ 


U\Dlf + \D2f) - - U\H+\" + - l\Ho\" - ^,Re{S^Bt) 

V 

9 

PlY§ 

[Re{H2D2* +H2D2*) + 2^/iRe{H2Dt)\ 

V 

10 

p2y§ 

[|P^P + |P^P+3|Po^P] 

V 

11 

Pl^Re{Yi) 

3 

^10 

^ Pe((P^ + H^)S^*) - ^ Pe((Pi + H2)Bt) + ^ Re{{B^ + B^)Ht) 

I 

V 

12 

Pl^Re{Yi) 

_ 3 
5V2 

Re{{H2 + H^)Ht) + y^l Re{{B% + Bt)S^*) + ^ Re{{B2 + Bt)Bt) 


V 

13 

Pl^Re{Yi) 

[2Re({Dl + Dt)Ht) + V3Re{{H2 + H2)Dt)] 

V 

14 

Pl^Re{Yi) 

-^Re{{D^+D^)Dt) 

V 

15 

Pl^Im{Yi) 

3 

V2 

^ Im{{H^ - H^)S^*) + 1 Im{{B2 - B^)Ht) - ^ ImHH^ - H2)Bt) 


V 

16 

P^V2Im{Y^) 

3 

^ ImiiBi - B2)B2*) + i ImiiH^ - Hb)Ht) + ^ Ira{{Bl - B2)S^*) 


V 

17 

PlV2Im{Yi) 

^ [2Im{{D2 - D2)Ht) + V3Im{{H2 - H2)D2*)] 

V 

18 

PlV2Im{Y^) 

1 ImiiD^ - D2)D2*) 

V 

19 

PSV2Re{Yi) 

+ ReiD^Dt*)] 

V 

20 

P2V2Re{Yi) 

-| [Re{H2D2*)+ Re{D2H2*)] 

V 

21 

P2V2Re{Yi) 

^ [-i Pe(P^P£*) + 1 ReiH^H^*)] 

V 

22 

P2V2Re{Yi) 

I yf [ReiH^Dt*) + ReiD^H^*)] 

V 

23 

P2V2Re{Yi) 

4 / 3 

TVs ^e(U+lJ_ ) 

V 

24 

P§V2Im{Yi) 

sjl [imiH^Ht*) + ImiDlD’l*)] 

V 

25 

P?V2Im{Yi) 

1 ImiH^Db* + DIh2*) 

V 

26 

P2V2Im{Yi) 

^ [i Im{DXD2*) - 1 Im{H2Hb*)] 

V 

27 

P2V2Im{Yi) 


V 

28 

P2V2Im{Yi) 

4 / 3 T,y^ { TiL ■r\L*\ 

“7 V 5 ) 

V 

29 

p2y° 

^ [{\h^\^-\h2;\^) + {\d^\^-\d2\^)] 

- (L^R) 

30 

p2y2 

Re{H^D’P - P+P+*) 

V 

31 

p2y2 

[f (|P^P - \Dl\^) - - IP^P)] 

V 

32 

p2y2 

- h2d2*) 

V 

33 

p2y2 

_2^(|^L|2 _ |^n|2^ 

V 

34 

PlV2Re{Y2) 

V 

^ [V2>Re{{Hl - H2)S^*) + ViRe({B2 - B2)Ht) - Re{{H2 - H2)Bt)] 

V 

35 

PiV2Re{Y2) 

3 

V2 

Re{{Hi - H^)Ht) + ^ Re{{Bl - B^_)S^*) + Ay^Pe((Pi - B^^)B2*) 


V 

36 

PiV2Re{Y2) 

^ [2Re((B^ - + V3Re((II^ - H2)Bt)] 

V 

37 

PiV2Re{Y2) 

IV5Re((B^-B^)Bt) 

V 

38 

P/V2/m(y/) 

-\/w lV5Jm((Ri + R2)S^*) + V3Im((Bi + B2)H2*) - Im{(Hl + Hh)B2*)\ 

V 

39 

PiV2Im{Y2) 

-\Jh W^Im{{H2 + H2)Ht) + V5Im{{B2 + B2)S^*) + f Im{{B2 + B2)Bt)] 

V 

40 

PiV2Im{Y2) 

-7h l2Jm((B2 + B2)II2*) + V3Im((I{2 + II2)B2*)] 

V 

41 

PiV2Im{Y2) 

-IV5Im((B2+B2)Bt) 

V 


TABLE III: The helicity-basis moments of the 41 orthonormal angular functions /i(n) in Eq. 46. 
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i 

m) 



1 

PSYo° 

[\H2\^ + |77|f p + |i7f p + + |Di|2 + |75i|2 

+ {L R) 

2 

p?yS 

2 -YRe{H2D2*) + Re{S^H2*) + Re{H^Df^* + T>i*) 



3 

p0\^0 

^ (IHf p + lllip) - (|7ff |2 + |77f |2) + \h^\2 + \d^\2 + 2 Re{S^D2*) 


4 

pO-^0 

-^3 ^0 

[- 7?e(77|f Df* + HiDi*) + ^7?e(77o^D^*)] 

V 

5 

pO’^O 
-^4 ^0 

f [-2(|7?fp + |Ili|2) + 3|Do^|2] 

« 

6 

PoY^ 

[(|7?f P + \Di\^) + (|77,[p + IHfp) - 2|S^|2 - 2|7?o^p - 2|J7o'^p] 

V 

7 


2^ 7?e(J7|f Df* + HiDi*) - ^ Re{S^H2*) - | ReiH^D^*) 


V 

8 

P^Y^ 

Y{\D^\^ + \Di\^) - I\D2\^ - - l\Ho\" " ^-.ReiS^D^*) 

« 

9 

pM 

[Re{H^D^*+H2D2*) + 2VSRe{H2D2n] 


10 

p2y^ 

[\D^\^ + \D2\^ + 3\D2\^] 

V 

11 

Pl^Re{Y2) 

3 

yio 

y'l 7?e(77|f S^*) - y^A 7?e(77f 77o^*) + y'f 7?e(7)f 77o^*) 


V 

12 

PiV2Re{Yi) 

3 

5 

7?e(77|f7fo'^*) + y^7?e(77^V^*) + ^ 7?e(77f77f*) 


V 

13 

PiV2Re{Yi) 

-SVTI [2Re{D^H2*) + V3Re{H;;D2n] 

« 

14 

PiV2Re{Yi) 



15 

PlV2Im{Y2^) 

3 

^ 7m(77f S^*) + 1 ImiDiH^n - ^ Im{HtD2*) 


« 

16 

PiV2Im{Yi) 

3 

^ Im{DiD2*) + 1 7m(77f 77o^*) + ^ 7m(77iS^*) 



17 

PiV2Im{Yi) 

^ [2Im{DiH2*) + V3Im{H2D2*)] 

V 

18 

PlV2Im{Yi) 

^ImiDiD^n 

V 

19 

P2V2Re{Yi) 

-2VIS [d^ii 1' - + (l^f P - \d2 \")] 

V 

20 

P2V2Re{Yi) 

-| [Re{H^D^*) - Re{D2H2*)] 

« 

21 

P2V2Re{Yi) 

# [-H\Di\^ - \Di\^) + P - \h2\^)] 


22 

PSV2Re{Yi) 

1^1 [i?e(J7fD[*) - Re{DiH2*)] 

V 

23 

P!V2Re{Yi) 


« 

24 

PSV2Im{Yi) 

y^l [7m(77i:H|f*)+ ImpiUf*)] 

V 

25 

P2V2Im{Yi) 

|/m(77fD|f* + Dii7|f*) 

V 

26 

P2V2Im{Yi) 

^ [i Im(D2D^*) - 2 Im{H2H^*)] 

V 

27 

P2V2Im{Yi) 

- 1 y^f + iLfDf) 

V 

28 

P2V2Im{Yi) 

4 Is r.™ / j~^L 7~jL*\ 

“7 Vs-' '-^2-■'^11 2 

V 

29 

PoY2 

-V3 [i?e(i7i;i7f*) + Re{D2D^*)] 

- (L^R) 

30 

p2y2 

-^Re{H2D^* + H^D2*) 

V 

31 

p2y2 

-Tm [f R<DiDf*) - Re{H2H^*)] 

V 

32 

p2y2 

^ Re{H2Dl(* + H^Di*) 

« 

33 

p2y2 

ReiDiD^*) 


34 

Pl^Re{Y2) 

Re{H2S^*) + V3Re{D2H2*) - ReiH^D^*)] 

V 

35 

P^^Re{Y2) 

3 

^ Re{HlH2*) + ^ Re{DiS^*) + Y^^Re(DiD2*) 


V 

36 

PiV2Re{Y2) 

^ [2Re{DiH2*) + V3Re{H2D2*)] 

V 

37 

P2V2Re{Y2) 

Y(p-Re{DiD2*) 

V 

38 

P2V2Im{Y2) 

-\ll + ^Im{D[^H2*) - /m(77|f D^*)] 

« 

39 

PiV2Im{Y2) 

-y^l [Vilm{HlH2*) + V^ImiDp^*) + f 7m(77[79f*)] 


40 

PiV2Im{Y2) 

-6^^ [2 7m(77f77o^*) + V3 7m(77|f7)^*)] 


41 

PiV2Im{Y2) 

-fVl0 7m(77f77^*) 



TABLE IV: The transversity-basis moments of the 41 orthonormal angular functions in Eq. 46. 
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from which, the efficiency-corrected true moments can be 
calculated as 


h = (54) 

Likewise, the measured covariance matrix of the moments 
is estimated as 


a. 




Eijkbk, (55) 


and the covariance matrix of the acceptance corrected 
moments are 


= {E-^UCki(,E-^)3i. (56) 

In the next step of the method of moments (MOM), 
if the moments functions are parameterized by a set of 
parameters a in some physics-motivated model as bi(d), 
the values of the a can be obtained by minimizing the 
function 

[dj - bjid)] (57) 

ij 


2. Unbinned maximum-likelihood fits (UML) 


In the equivalent unbinned maximum-likelihood 
(UML) method, the efficiency incorporated pdf is 


g'{VL,d) = 


Y,h{d)Mn)ein) 

^bi{d)Ei 


(58) 


The likelihood function to maximize is 

^ a.t 


E{d)= g'{VLk,d), 


k^l 


(59) 


which leads to the negative log-likelihood (NLL) to min¬ 
imize as 


-21n(/:(d)) = 


'^b^{d)Ei 


AT" 


-2 In 


k^l 


^h{d)f^{nk) 


(60) 


B. Studies with toy Monte Carlo 


To validate the above expressions, we consider a simple 
rate expression for toy studies: 

dN N 

-—= g{0, a, 13) =a cos 9P sin 9) (61) 

d9 TT -\- 2 p 


Set 

efficiency t{6) 

I 

(1 -tsin26»)/2 

II 

(1 ± cos“ 9)/2 

III 

(1 -tcos^6l)/2; 0 for |6» - f 1 < 0.1 


TABLE V: The different efficiency functions used in the toy 
studies of Sec. IX B. 


with 9 G [0,7r] and {a, (3} being the target parameters 
to be determined. The total number of events, N, is a 
nuisance parameter for the moment. The orthonormal 
basis functions are 


/i = ^ 

Vtt 

(62a) 

COS0 

vT2 

(62b) 

sin0 —2/7r 
\/7r/2-4/7r’ 

(62c) 

and the corresponding moments 


bi — 

Vtt 

(63a) 

Na^JVi2 

^±2/? 

(63b) 


(63c) 


Without any loss of generality, we model the detector 
efficiency as the three sets of functions given in Table V. 
Specifically, we note that Set III incorporates a “hole” in 
the detector around 0 = 7r/2, where the efficiency drops 
to zero. 

Figure 5 shows the pull distributions from fits to a 
thousand toy samples generated according to Eq. 61 and 
the efficiency function as Set III in Table V. 


C. Incorporating background 

Next, to incorporate background, we assume that there 
is a discriminating variable m, un-correlated with the 
angular variables O. Let there be events (signal 

and background combined) in the “signal region” in the 
variable to, and N'° events in a suitably defined “sideband 
region”, containing pure background events. Also, let hb 
be the estimated background under the signal peak in 
the “signal region”, obtained from a signal-background 
separation fit in the variable to. 

Independent toy sample sets with different purity levels 
were generated. Figure 6 shows the case for a toy sample 
with the discriminating variable to = mb representative 
of the B mass. The signal lineshape is a Gaussian while 
the background is constant. The “signal region” is chosen 
as ±2(7 around the mean, as obtained from the signal- 
background fit. The low and high sideband regions are 
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pull In a, MOM 




-4 -2 0 2 4 

pull in |3, MOM 



pull In a, UML 


pull In (3, UML 


FIG. 5: (Color online) Pull distributions from a thousand toy 
samples of fits to Eq. 61 with the efficiency function as Set III 
in Table V. The upper and lower plots use the MOM and 
UML techniques, respectively. No background is included. 



uIb (MeV) 


FIG. 6: (Color online) A toy sample incorporating both ef¬ 
ficiency and background effects used in the validation study. 
A fit to the profile is shown as well. 


negative weights to the events in the sideband region: 


-2ln{C'{d)) 


fib) In 


biid)E^ 


^data 

-2 y] In 

k^l 

+ 2 x In 


J2b^{d)M^k) 

. i 

J2b^{d)Mnk) , 


(64) 


where x = fi^/N-b is a scale factor relating the back¬ 
ground level under the signal, to that in the side-band. 

Following the derivation in Refs. [11, 23], the covari¬ 
ance matrix from minimizing the pseudo-likelihood func¬ 
tion in Eq. 64 has to be modified to yield the true covari¬ 
ance matrix, C'*’, incorporating the additional uncertain¬ 
ties due to the background subtraction part as: 

= C[l + {n^{l + x)g + NlalH]C], (65) 

where 


Wb 


IVb ^ 




51ng'(flfc,a) 91ng'(rjfe,a) 


dax 


da. 


( 66 ) 


1 a In g'(fife, a) 
iVb dax 


1 d\ng'{ni,d) 

Nh ^ da^j, 

(67) 


and C is the covariance matrix returned by the HESSE 
routine. Summing over repeated indices, the partial 
derivatives are explicitly 


d\ng'{nk,d) 

dax 


d b,{d) r /i(flfe) _ Ej 

dax [bj(d)fj{nk) bj{d)Ej_ 

( 68 ) 


and Gx is the uncertainty on the background scale factor 

X. 

In the moments expansion method, the background- 
subtracted measured moments and the covariance matrix 
are estimated as 

jydata jyb 

M^k)-xJ2M^k) (69) 

k^l k^l 

^^data 

= E M^k)fj{nk) + E M^k)fj{^k) (70) 

k=l k=l 


The pull distributions from the MOM and ULM fits 
and the corresponding covariance matrices and C^, 
respectively, are shown in Fig. 7. 


taken as tob < 5200 MeV and ms > 5360 MeV, respec¬ 
tively. The background is generated flat in and 0, but 
folded with the relevant efficiency functions in Table V. 

The “pseudo-likelihood” is then defined by assigning 


D. Discussion 

We point here to some of the salient features of the 
MOM. The set of moments in Eq. 46 constitute a concise 
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written in terms of the moments as: 


n = 


(ri + v^r3)(ri-v^r3/2) 


35- 


(71) 


The important point to note here is that no complicated 
multi-dimensional angular fit is required for any of these 
checks. 

We would also like to comment on the use of the nor¬ 
malization integrals in Eq. 52 as opposed to analytic 
modeling of the efficiency function and reweighting of 
events by the inverse of the efficiency. The latter involves 
a complicated fit which can be unstable without due to 
local “holes” in the acceptance function. The normaliza¬ 
tion integrals, on the other hand, are found to be more 
robust under these situations. 


FIG. 7: (Color online) Pull distributions from a thousand 
toy samples of fits to Eq. 61 with the efficiency function as 
Set III in Table V. The upper and lower plots use the MOM 
and UML techniques, respectively. The samples incorporate 
a signal to background ratio corresponding to that in Fig. 6. 


representation of all the angular information content in 
the entire dataset. The relations between the different 
moments and the amplitudes are ab initio not built in. 
These relations can be used as checks for understanding 
of the detector acceptance. They can also be incorpo¬ 
rated during the model-dependent minimization fit as 
described by Eq. 57. If the model-dependence is reliably 
known, the MOM and a direct UML fit give the same 
results, as we explicitly demostrated in Sec. IX B. 

However, if the underlying physics model is unknown, 
the MOM can provide simple and model-independent 
confirmations of certain interesting physics features. For 
example, as pointed out in the introduction, a complex 
RH admixture cr in the weak hadronic current leads to 
angular terms proportional to siny in SL decays, that 
are absent in the SM. The presence of these terms in the 
data can be examined using any of the moments in Ta¬ 
ble IV corresponding to Im(Y^), where m 0. If the 
statistical significance of these moments are found to be 
high enough, this could constitute tension with the SM. 

Similarly, the observables (|i7oP + |<S'P), |77 {||._l}P, 
l^{o.||,-L}P individually expressed in terms of the 

moments in Table IV. Therefore, if one is interested in 
the presence of a D-wa.ve component under the K*{892) 
for H —)• RT this can be directly probed via the 

moments. In the absence of a D-wave component, the 
observables |i7oP and |S'p can also be extracted directly 
from the moments, allowing an estimate of the S'-wave 
fraction. For the observable Pg [25] that is predicted to 
be theoretically clean at low the LHCb collaboration 
has recently observed [26] a 3.7a deviation from the SM. 
In the absence of non-P-wave components, this can be 


X. SUMMARY 

In summary, we provide expressions for the full angular 
decay rate in P —X.^i .^2 decays where the £2 lepton can 
be either a charged {e,/r} or a neutrino. We considered 
the final state X to include complex S-, P, and P-wave 
amplitudes. The rate expression is expanded in a basis 
of orthonormal moments functions and a procedure to 
extract the corresponding moments employing a counting 
measurement is desribed and validated. We expect the 
present work to be directly applicable to ongoing analyses 
at BaBAR and LHCb. 


Appendix: Angle definitions 


In this appendix we provide the explicit definition of 
the angles in terms of the 3-vectors. The definitions are 
equivalent to the GS definitions as explained in Sec. HB. 

We follow the convention adopted in App. Ref. [13] 
that the superscript on any 3-vector denotes the ref¬ 
erence frame. For any ordered four-body final state 
B -A {Pi,P 2 ,^ 1 ,^ 2 } where P{i, 2 } are pseudoscalars and 
^{ 1 , 2 } are leptons, we define 


P^1^2 = Pii + W2 
Qiit2 = Ph - Pi2 
PPlP 2 = Pil +PP 2 
QPiP2 = Ptl -PP2- 


The helicity angles are defined as 


cos 61 


cos 6 v 


Ql 


•p; 


P1P2 


\Qi[ 




(jPP 


\pU I 

PPlP2\ 

pPP 


iQ?rP2i 


ppp I 


(A.la) 
(A.lb) 
(A.lc) 
(A.ld) 


(A.2a) 


(A.2b) 
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where U and PP in the superscripts refer to the leptonic 
and hadronic rest frames. 

The normals to the two planes are defined as 

^ Qfii 2 (A.3a) 

lVpiP2 = Pp^p^ X Qp^p^, (A.3b) 

and the dihedral angle between the planes is defined by 


leading to a single sign flip in cos0i compared to the 
EWP theory convention, as was explained in Eq. 9. 

For the CP conjugate decay B —> {A.P 2 ^i^ 2 }, we per¬ 
form the CP conjugation explicitly while maintaining the 
order. The same procedure using Eqs. A.l- A.4 is applied 
to the CP conjugated system to yield the angles. This 
leads to a single sign flip in the angle Xi mentioned 
earlier in Sec. VI. 


cosx = 


smx = 


Neie2 ■ A^PiP2 
l^^l<!2ll^PlP2l 

X ^PlP2 \ 

j^c^2ii^PiP2i;' \pF,iS 


p. 




(A.4a) 

(A.4b) 


For the B decay, our ordering is 5° —> {K 7r+^ 1'+}, 
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